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Abstract 



We investigate the spatial overlap of nonclassical ultrashort light pulses 
produced by self-phase modulation effect in electronic Kerr media and its 
relevance in the formation of polarization-squeezed states of light. The light 
polarization is treated in terms of four quantum Stokes parameters whose 
spectra of quantum fluctuations are investigated. We show that the frequency 
at which the suppression of quantum fluctuations of Stokes parameters is the 
greatest can be controlled by adjusting the linear phase difference between 
pulses. By varying the intensity of one pulse one can suppress effectively the 
quantum fluctuations of Stokes parameters. We study the overlap of nonclas- 
sical pulses inside of an anisotropic electronic Kerr medium and we show that 
the cross-phase modulation effect can be employed to control the polarization- 
squeezed state of light. Moreover, we establish that the change of the intensity 
or of the nonlinear phase shift per photon for one pulse controls effectively 
the squeezing of Stokes parameters. The spatial overlap of a coherent pulse 
field with an interference pulse produced by mixing two quadrature-squeezed 
pulses on a beam splitter is analyzed. It is found that squeezing is produced 
in three of the four Stokes parameters, with the squeezing in the first two 
being simultaneous. 
PACS: 42.50.Dv, 42.50.Lc. 

Keywords: ultrashort light pulses, electronic Kerr nonlinearity, self-phase modulation, 
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I. INTRODUCTION 



Quantum states of light with squeezed polarization have recently been in the focus of 
theoretical and experimental investigations. The quantum analysis of the light polariza- 
tion is generally based on four Stokes parameters associated with four Hermitian Stokes 
operators [1,2]. By definition, the quantum state of light with the level of quantum fluc- 
tuations of Stokes operators smaller than the level corresponding to the coherent state is 
called polarization-squeezed (PS) state. Its existence was predicted theoretically in [3] and 
was revealed experimentally in a series of recent experiments [4-6]. Historically, the first 
experimental realization of a PS state of light implied the mixing of a strong orthogonally- 
polarized coherent beam with squeezed vacuum on a 50/50 beam splitter (BS) [4]. The 
optical parametric amplification (x^) has also been successfully employed to generate the 
PS state of light [7-9]. In the last decade the improvement of fiber optics techniques in the 
pulse field regime facilitated experiments directed to produce nonclassical solitons [10,11] 
or highly stable zero-dispersion quadrature-squeezed ultrashort fight pulses (USPs) [12-14]. 
Recently it has been found experimentally that the simple spatial overlap of two orthogonally 
polarized quadrature-squeezed pulses leads to the formation of the PS light [5,6]. 

The self- phase modulation (SPM) effect, which is responsible for the quadrature squeez- 
ing of pulses in optical Keer-fike media (x^), does not infiuence the photon statistics of 
pulses since the photon number is a constant of motion [15]. It was noted for the first time 
in [16], that a consistent time-dependent quantum description of the SPM effect must nec- 
essarily account for the additional noise related to nonlinear absorbtion (the imaginary part 
of nonlinear susceptibility). In [17] the Kerr nonlinearity has been treated as a Raman- like 
one, the quantum and thermal noises being considered as a fiuctuating addition to the re- 
laxation nonlinearity in the interaction Hamiltonian. The addition is necessary in order to 
preserve the commutation relations among the pulse field-amplitude operators. Hence, the 
resulting quantum equation of motion for the total field includes the Hermit ian-phase noise 
operator which describes the coupling of the field to a collection of localized, independent, 
medium oscillators. Besides, the average value of the noise operator on the coherent state 
is considered to be zero. The time-delayed Raman response of nonlinearity, which is around 
50 fs in fused-silica fibers, varies over frequencies of interest and is caused by back action 
of nonlinear nuclear vibrations on electronic ones. However, the contribution of Raman 
oscillators to the Kerr effect is secondary to the one of the electronic motion. Indeed, if we 
deal with the nonlinear propagation of USPs, for instance through fused-silica fibers, since 
the Raman oscillators contribute with less than 20% to the Kerr effect, the contribution of 
the electronic motion on ~ 1 fs time scale is more than 80% [18]. The analysis of the SPM 
effect then has to focus on contribution of the electronic motion, especially in the case of 
USPs with a duration much less than the time-delayed Raman response or in case the pulse's 
frequencies are far from any Raman resonances. An attempt to develop a quantum theory 
of the pulse SPM primarily due to the electronic motion has been undertaken in [19], but 
the electronic nonlinearity has been modeled as a Raman-like one. The electronic nonlinear 
response function in the frequency domain over the bandwidths of interest is considered 
constant and approximated with a delta function in the time domain. This implies an in- 
stantaneous electronic response of the Kerr nonlinearity. In real situations the nonlinear 
response of electronic nonlinearity is finite. It results in a non-delta relaxation function in 
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the proper analytical description of the Kerr effect. Such a consistent description based on 
the momentum operator, which is connected with the evolution of the field in space and 
incorporates the relaxation function of the electronic nonlinearity, has been implemented 
in [20,21], where the correlation functions and spectra of quadrature components of USPs 
subjected to the SPM effect in the electronic Kerr medium are investigated. Since the 
momentum operator is in the normally ordered form, in order to satisfy the commutation 
relations for time-dependent Bose operators there is no need for additional fluctuating terms. 
The approach developed in [20,21] for the SPM effect has been extended in [22] to the case 
of two USPs undergoing, besides the SPM effect, the cross-phase modulation (XPM) effect. 

In this paper we apply the quantum theory of the SPM and combined SPM-XPM effects 
developed in [20,21] and [22], respectively, to the case of spatial overlap and interference on 
a BS of quadrature-squeezed pulses. The nonclassical USPs are obtained by undergoing the 
SPM effect in an inertial electronic Kerr medium. The finite relaxation time of the electronic 
Kerr nonlinearity is accounted for and the dispersion of linear properties is described in the 
first approximation of the dispersion theory. We show that the approach for the SPM effect 
based on the momentum operator allows an adequate calculation of correlation functions 
and fluctuation spectra of quantum Stokes operators. 

The paper is organized as follows. In Section II the quantum model of SPM effect for 
USPs based on the momentum operator for the pulse field is exposed briefly. We present 
the space-evolution equation for time dependent Bose operators, its solution, as well as 
some elements of the algebra of Bose operators. Section III introduces the quantum char- 
acterization of polarization of light in terms of Stokes operators, as well as their correlation 
functions and spectra. In Section IV we analyze the overlap between a coherent USP and 
a quadrature-squeezed one, the latter USP being produced by employing the SPM effect in 
an electronic isotropic Kerr medium. The average values, the correlation functions of Stokes 
operators ^2 and 5*3 are computed, and the spectrum of 5*2 is graphically represented. This 
analysis is continued in Section V by investigating the spatial overlap of two independent 
quadrature-squeezed USPs. In Section VI we study the overlap of two quadrature-squeezed 
pulses inside a nonlinear anisotropic Kerr medium in the presence of the combined SPM- 
XPM effect. The momentum operator for the XPM effect introduces the section. Here we 
show that the XPM effect can be employed to control the formation of polarization-squeezed 
spectra. Section VII is dedicated to the interference on a BS of two independent quadrature- 
squeezed USPs and their overlap with a coherent USP or with squeezed vacuum. We reveal 
that the squeezing is possible in three Stokes parameters 5*0, Si, and 5*2. We compare our 
analytical results with the experimental ones reported in [5,6]. Our concluding remarks are 
presented in Section VIII. 



II. SPM EFFECT IN ELECTRONIC KERR MEDIUM 

The traditional way to derive the quantum equation for the SPM effect is based on the 
interaction Hamiltonian. As a consequence one obtains a time-evolution equation. The 
transition to space-evolution equation is usually realized by enforcing the replacement t 
z/u, where z is the distance passed inside the Kerr medium and u is the group velocity. This 
approach is reasonable for single-mode radiation. However, the analytical description of the 
nonlinear pulse propagation contains both t and z. The consistent quantum description of 
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the Kerr effect requires tlien tlie use of the momentum operator connected with the evolution 
of the pulse field in space [23] . 

When accounting for the finite relaxation time of the Kerr nonlinearity the SPM effect 
is described with the following momentum operator [20-22]: 

^ roo rt ^ 

Gspm(z) ^hp dt H(t - ti)'N[h(t, z)n(ti, z)]dti (1) 

J — OD J—OO 

where h is the Plank constant, the factor (3 is defined by the Kerr nonlinearity of the medium 
{P — fiwokon2/SnoV [22,24]), N is the normal ordering operator, H{t) is the causal nonlinear 
response function [H{t) 7^ at i > and H{t) = at t < 0], n{t,z) = A+{t, z)A{t, z) is 
the "photon number density" operator in the cross-section z of the medium, and A(t, z) 
and z) are the photon annihilation and creation Bose operators with the commutation 

relation [A(ti, ^), A+(t2, ^)] = 5(^2 — ^1), respectively. The thermal noise is neglected and 
the expression (1) is averaged over the thermal fluctuations. In our approach the pulse 
duration Tp is much greater than the relaxation time r,., and the Kerr medium is lossless and 
dispersionless; that is, the frequencies of the pulse are off resonance. 

In Heisenberg representation, the space evolution of A{t, z) is given by the equation [23] 



. JA{t,z) 

a 

oz 



G,p^iz),A{t,z) . (2) 



Then, with (1) the space evolution equation of A{t,z) in the moving coordinate system 
{z — z , t — t — z/u, where t is the running time and u is the group velocity) 



has the solution given by 



^^-dit)Ait,z) = (3) 



i(t,^) = e«Wi(t), (4) 



where 0{t) = i^q{n{t)), 7 = j3z, h{t) = A~^{t)A[t) is the "photon number density" operator 
at the entrance into the medium, n(t) = n{t,z = 0), and q{n{t)) = J!^^h{ti)n{t — ti)dti, 
with h{t) = H{\t\). 

In comparison with the so-called nonlinear Schrodinger equation, used in the quantum 
theory of optical solitons (see for instance [25] and references therein), in (3) the light 
pulse dispersion spreading in the medium is absent. This approach corresponds to the first 
approximation of the dispersion theory [24]. Note that the structure of q{n{t)) is like that 
of the linear response in the quantum description of a USP spreading in the second-order 
approximation. 

The statistical features of pulses at the output of the medium can be evaluated by using 
the algebra of time-dependent Bose operators [20,21]. For example, in such algebra we have 
a permutation relation of the form: 

i(i,)e^(*2) = e^(*^)+^(*^-*i)i(ii) (5) 

where 'D{t2 — h) — i'yh{t2 — ti), 27 is the nonlinear phase shift per photon. By using the 
theorem of normal ordering [20,21] one obtains the average values of the Bose operators over 
the initial coherent summary state \cy{t)): 
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^gO(t)^ (6) 

The parameters (/){{) = 2^fn{t), fi{t) = j'^n(t)/2 are connected with SPM of the pulse. Here 
(f){t) is the nonhnear phase addition caused by SPM. The eigenvalue a{t) of the annihilation 
operator A{t) over the coherent state \a{t)) can be written as a{t) — \a{t)\e^'^^*\ where (p{t) 
is the linear phase of the pulse. Then |Q;(t)p = {n{t)) = n{t). The time dependence can be 
separated in \a{t)\ by introducing the pulse's envelope r(t) so that \a{t)\ = \a{0)\r{t) with 
r(0) = 1. For the simplicity of notations let 0(t = 0) = 0o and n(t = 0) = no- In (6) and 
(7) /C(ti,t2) = IJ'O^'^iii + ''"/2)5'(r) is the temporal correlator due to SPM, where iiq = ■yno, 
g{T) = (1 + |t|/t^)/i(t), and r = ^2 - h. 

Notice that in our approach there is no assumption about the form of the relaxation 
function H{t) of the Kerr nonlinearity. It is positive at any instant t > and zero otherwise. 
If the SPM efect is due primarily to the electronic motion occurring on ~ 1 fs time scale 
{jr < 1 fs), then in the absence of one- and two-photon and Raman resonances, the relaxation 
function at i > can be approximated in the form [24]: 

H{t)^{l/Tr)e^p{-t/Tr). (8) 

In the following sections this specific form of the relaxation function is used when evaluating 
the Fourier transform of /i(t) and g{T). 



III. QUANTUM CHARACTERIZATION OF THE LIGHT POLARIZATION; 

STOKES OPERATORS 

The quantum description of the polarization of light is usually done in terms of four 
Stokes operators Sa {a — 0, 3) [1,2]. For two spatially overlapping USPs the Stokes operators 
are introduced as follows: 

Soit) = Atit)A^{t) + i+(t)i2(t) = Tr[Dtit)D^it)], (9) 

S,{t) = Atit)A^(t) - Atm^it) = Tr[a^b+{t)Di{t)l (10) 

S^it) = i+(t)ii(t) + At{t)A^{t) = Tr[a,bt{t)b,{t)i (11) 

S^{t) = i[A+{t)A,{t) - At{t)A2{t)] = Tr[-a2b+{t)b,{t)] (12) 

where 0"^ are Pauli matrices (A; = 1, 3), and I?-matrices are: 

\ A,{t))' \A,(t) ) ^ ' 

The definition of Stokes operators Sk in terms of Pauli matrices (10)-(12) suggests that in 
quantum optics they play a similar role to the Pauli matrices dk which in quantum theory 
of angular momentum describe the rotations in two-spinor formalism. Indeed, apart from a 
normalization factor, the commutation relations for Stokes operators are identical to those 
for Pauli matrices; [Si{t), Sj{t)] = 2ieijkSk{t), where indices k = 1, 2, 3 are taken by cyclic 
permutations. The energy operator S'o(t) commutes with any Sk{t), i.e., [So{t), Sk{t)] — 0. 
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In addition, the Stokes operators generate a special non-abelian unitary group of symmetry 
transformations SU2 that obeys the same Lie algebra as the three-dimensional rotation group 
of Pauli matrices. Since the Stokes operator 5*0 (t) describes the total energy of the two pulse 
field at the time t, Sk{t) characterize light polarization and form a Cartesian axis system. 
The polarization state is visualized as a vector on the Poincare sphere and each point on 
the sphere corresponds to a definite polarization state whose variation is characterized by 
the motion of the point on the sphere. If the Stokes vector points in the direction of Si, S2, 
or 5*3, than the polarized part of the beam is horizontally, linearly at 45°, or right-circularly 
polarized, respectively. When Stokes vectors of two USPs point in opposite directions pulses 
do not interfere. The radius of Poincare sphere, i.e., the average length of the Poincare 
vector 7^(t) = (Xl=i{Sk{t))^y/^, defines the average intensity of the polarized part of the 
radiation. The ratio of the intensity of the polarized part to the total average intensity 
{So{t)) is called degree of polarization V{t) = ^.{t) / {So{t)) and it is an important measure 
in quantum optics. 

To analyze the behavior of quantum fluctuations of Stokes operators we introduce 
their correlation functions: 

Rs^{t,t + T) = {Sa{t)Sa{t + T)) - {Sa{t)){Sa{t + T)). (14) 

In the frame of the quantum model for the SPM effect presented in Section II the correlation 
functions (14) can be analytically computed by using permutation relations of the type 
(5) and the averages (6) and (7) of time-dependent Bose operators (see for example [26]). 
The average values in (14) are calculated on the summary coherent quantum state \a{t)) = 
\ai{t))®\a2{t)) . We assume here that the initial USPs are in coherent states, the j-th pulse's 
operator Aj{t) acting only on the corresponding state vector \aj{t)) within the associated 
sub-Hilbert space Tij, so the factorization of the quantum states of the two sub-Hilbert 
spaces takes place. The orthogonal basis of summary states belongs to the global Hilbert 
space Ti = Tii ® 7^2- This basic assumption is followed in the paper when calculating the 
quantum average of physical quantities of interest. If both USPs are in coherent states, 
then Rs^{t,t + r) = 5(r). As is well known from experiments, if USPs undergo some 
nonlinear transformations (for example, parametric ampliflcation or SPM effect) and further 
overlap spatially, than the formation of PS state of light is allowed. This means that, in the 
theoretical description of the formation of the PS state in such processes, the correlation 
functions of Stokes parameters are different from 5-function. However, the polarization 
squeezing cannot be obtained in all Stokes parameters at the same time. We show bellow, by 
analytically deriving the correlation functions of Stokes parameters that the spatial overlap 
of quadrature-squeezed pulses, the overlap of nonclassical pulses inside an anisotropic Kerr 
medium in the presence of combined SPM-XPM effect, or the interference of nonclassical 
USPs on a BS, lead to the formation of the PS state. 

With the correlation function calculated in agreement with (14) one can proceed with 
the evaluation of the spectra Ss„{uj,t) of quantum fluctuations of Stokes operators at any 
instant t by merely applying the Wiener-Khintchine (WK) theorem, which states that 

SsA^,t) = r RsSt.t + T)e'^^dT. (15) 
J —00 

For two coherent USPs that overlap spatially the spectrum of quantum fluctuations of 
is frequency independent and constant, i.e., Sg'^{u!,t) — 1. Consequently, the deviation 
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from the coherent level Sg^{u!,t) — Ss„{u!,t) — 1 can be used as a measure of quantum 
fluctuations behavior. Thus, since the case < Sg^^{uj,t) corresponds to anti-squeezing, the 
case — 1 < Sg^{u!,t) < corresponds to the suppression of quantum fluctuations of Sa- 

IV. OVERLAPPING COHERENT AND QUADRATURE-SQUEEZED PULSES 

In the quantum model for the SPM effect in Kerr media developed in [20,21] and shortly 
exposed in Section II the finite relaxation time of the electronic nonlinearity is accounted. 
Note that this approach leads to an accurate calculation of correlation functions and fluctua- 
tion spectra of quadrature components for USPs subjected to the combined SPM-XPM effect 
in electronic Kerr media [22]. Within the quantum model [20,21] we consider the spatial 
overlap of a coherent USP {Ai{t)) with a quadrature-squeezed one {A2{t,l)). The nonclas- 
sical USP is obtained by passing a Kerr medium of length I [7 = pi, A2{t, I) — e'^(*M2(t)]. 
The Stokes operators are defined in agreement with (9)-(12) with D-matrices in the form: 

We point out that the Kerr effect does not affect the photon statistics [15] and, in the 
assumption that one can neglect the dispersion effects, the operators 5*0 and Si, as well as 
their dispersions, are conserved. Below, by investigating the fluctuation spectra, we show 
that the overlap of USPs produces the squeezing of quantum fluctuations in 5*2 or §3. With 
the use of Eqs. (4) and (6) wc estimate the average values of these operators. We have: 
{So{t)) = Eii ni{t), {Slit)) = niit) - Mt), 

(S^it)) = 2[ni{t)n2{t)]'/'e-^^^'^ cos [^^{t) - Vi{t)], (17) 
{Ss{t)) = 2[rii(i)ri2(i)]V2e-^^(*) sin [^^{t) - <Pi{t)l (18) 

where $2(^) = <p2{t) + ^2{t). Note that {S3{t)) is shifted in phase with — 7r/2 in comparison 
with {S2{t)), in complete agreement with the Heisenberg uncertainty relation for noncom- 
muting observables. To investigate the behavior of quantum fluctuation of Stokes operators 
iS'2 and iS'2 we proceed by evaluating their correlation functions (14). It is not difficult to show 
that in the approximation 7 ^ 1 by using Eqs. (4)-(7) the correlation function Rszit, t + r) 
has the following analytical form: 

Rs,{t,t + T) = 5{t) + h{T)niit)Mt) sin2[v9i(t) - <l>2(t)] 

+ g{T)ni{t)cj>l{t) sin2 [^i{t) - ^^{t)]. (19) 

With the correlation function (19) the fluctuation spectrum of §2 can be obtained by ap- 
plying the WK theorem (15). Allowing for a small change of the envelope r(t) during the 
relaxation time and accounting the fact that for the electronic nonlinearity (8) we have: 
JZo h{T)e''^^dT = 2L{lj) and J^^ g{T)e''^^dT = 4L'^{uj) [20-22], for the fluctuation spectrum 
of ^2 we finally obtain 

Ss,{n,t) = 1 + 2Lin)ni{t)Mt) sin2[(^i(t) - $2(t)] 

+ 4:L^{n)ni{t)cj)l{t) sin" [<pi{t) - $2(^)], (20) 
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where L{Q) = 1/(1 + fl"^). Here the dimensionless frequency fl = ujTj. is introdTiccd for the 
simpUcity of notations. The correlation function of S^, as weU as its fluctuation spectrum, 
can be easily obtained by shifting (19) and (20) in phase with tt/2, respectively. Since 
the spectral density of quantum fluctuations (20) depends on the relaxation time of the 
electronic nonlinearity, the latter deflnes the level of quantum fluctuations below to the one 
corresponding to the coherent state. The second term on the r.h.s. of Eq. (20), which can be 
negative for a definite linear phase difference A(p(t) = ^P2it) — y?i(t) between USPs, indicates 
clearly that the level of quantum fluctuations of 5*2 can be less than that corresponding to 
the coherent state, S'|P'*(n,t) ~ 1. Let A.ip{t) at a deflnite reduced frequency f2o has the 
form 

A^^,(t) = I arctan (^^^^^) - Mt)- (21) 

Then, at this Aipopt{t) the fluctuation spectrum (20) reaches the minimum value 

Ss,(no,t) = l + 2ni(i)0^(i)L2(Qo) -2rii(t)02(O^W[l + 02W^'(^^o)]'/'. (22) 
At any frequency fl the fluctuation spectrum writes: 

Ss,{Q,t) = 55.(f^o,t) + 2ni(t)02(t)[L2(f]) -L2(fio)] -2ni(t)02(t)[/>(5^) -i^(^^o)] 

x{l + (j)l{t)L{no)[L{n) + L(f]o)]}[l + 0^(t)L2(f])]-V2 (23) 

As already pointed out, in order to characterize the deviation from the coherent level we 
introduce the normalized spectral variance Sg,^{Q,t) = [Ss2i^,t) — l]/ni(t). In Figs. 1 and 
2 we display Sg^{Q,t) for different values of the nonlinear phase addition 02 (^)- For linear 
phase difference optimized at low frequencies f2o = the squeezing in 5*2 is maximal around 
Q (Fig. 1). In case Qo — 1, the squeezing take place essentially at high frequencies 
Q 1 {uj — 1/Tr) (Fig. 2). Thus, the choice of hnear phase difference allow us to produce 
the squeezing basically at the frequency of interest. Moreover, the increase of the nonlinear 
phase addition 02 (^) produces a better level of squeezing. 



V. OVERLAPPING QUADRATURE- SQUEEZED PULSES 

Now we focus our attention to the case of spatial overlap of two independent quadrature- 
squeezed USPs. The nonclassical pulses obtained at the exits of electronic Kerr media of 
lengths li and I2 further overlap spatially. The Stokes operators for the nonclassical USPs 
are defined in agreement with (9)-(12) in terms of annihilation [Aj(t,lj) = e*^^*^*Mj(t)] and 

creation [Aj'{t,lj) = Aj'{t)e'^i *^*^] operators (j = 1,2) with D-matrices in the form: 

Proceeding as in the previous section for the average value of 5*2 and 6*3 we have: 

{S2{t)) = 2[ni(t)n2(t)]'/'e-^^W-^^W cos [<|.2(t) - <I>i(t)], (25) 
{Ss{t)) = 2[ni(i)n2(i)]^/2e-'^^W-'^^« sin [^^{t) - $i(t)]. (26) 
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FIG. 2. As in Fig. 1 but for 00 = 1- 
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The correlation function of 5*2 calculated by using the relations (4)-(7) in the approximation 
Ij < 1 is 

Rs,{t,t + t) = 6{t) + /i(r)[ni(t)02(t) - Mt)^^)] sin2[$i(t) - $2(t)] 

+ ^(r)[ni(t)02(t) +n2(t)0?(t)]sin2 [<fi(t) - $2(t)]. (27) 

Compare the relations (17)-(19) and (25)-(27). Instead of we have now = 

4>i{t) + (pi{t) where 4>i{t) is the nonlinear phase addition due to SPM effect on the pulse 
1 after passing the Kerr medium of length h. With the WK theorem, for the fluctuation 
spectrum of 5*2 we get: 

Ss,{n,t) = l + 2L{n)[n,{t)Mt)-Mt)Mt)]sm2[^i{t)-Mt)] (28) 
+ AL^{n)[m{t)<Pl{t) + n2ml{t)] sin^ [$i(t) - ^,{t)]. (29) 

Besides adjusting the nonlinear phase differences (pjit), we can control the spectral density 
(29) by changing the intensities nj{t) of USPs. If the linear phase difference Aip{t) — 
ip2{t) — <Pi{t) at a definite frequency flo is taken in the form 

^^opt = t; arctan u\a2u\ , - u\a2u\^ + Mt) - W), 30) 

2 VL(Qo)[ni(i)(/)2(i) +n2(i)0i(i)]y 

then the fluctuation spectrum (29) reaches the minimum value 

Ss,{n^,t) - l + 2[ni(t)02(t) +n2(t)0?(t)]i:'(l^o) -2L(no){[rii(i)02(i) -fl2(i)</'i(i)]' 

+ L2(Oo)[ni(t)0^(t) +n2(i)0?(i)]^}^/^. (31) 

At any frequency VL the fluctuation spectrum is: 

SsMt) = Ss,{^o,t) + 2[ni(t)02(t) + n2{t)(t>\mL\^) - L\n^)] - 2[L(fi) - L{no)] 
x{[ni(i)02(i) - n2{t)ct>^{t)f + L{Slo)[L{Q) + L{no)][n,{t)4>l{t) + n2{t)ct>l{t)f] 

x{[n^{t)(t>2{t) - n2{t)Mt)? + L\no)[niml{t) + n^mUt)]''}''^' . (32) 

In Figs. 3 and 4 we present the normalized spectral variance Sg,^{Q, t) = [Sg'J^{fl, t) — l]/ni{t) 
for different relations between pulses intensities. One can see that, by increasing the intensity 
of the control USP 2, we achieve a better level of squeezing. 



VI. OVERLAPPING QUADRATURE-SQUEEZED USPS IN ANISOTROPIC 

KERR MEDIA 

When two USPs with orthogonal polarization and/or different frequencies overlap spa- 
tially inside an anisotropic Kerr medium, beside SPM effect of pulses, the XPM effect and 
parametric interaction occur. The parametric frequency conversion can be neglected in 
the assumption of a large phase mismatch |A| > min{/3|Ai(0)|, /3|A2(0)|} and Az ^ 1 
[3,24]. In this situation the XPM effect can be employed to control the spectra of quantum 
fluctuations of quadratures [22] , as well as those of quantum Stokes parameters [26] . The 
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FIG. 3. Normalized spectral variance Sg^{^, t) for linear phase difference Aip(t) chosen optimal 
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consistent description of the combined SPM-XPM effect is based on tlie momentum operator 



^ ^ fOO ft 

Gxpm{z) = hp dt H{t - t]) [hlit, z)h2{ti, z) + hi{ti, z)h2{t, z)] dti. (33) 

J—oo J—oo 

In the above expression j = (3z is the nonlinear coupling coefficient [22,26]. With the total 
momentum operator G{z) we derive the solution of the space-evolution equation (3) for the 
pulse 1; 

= e^iW+^2(*)ii(t), (34) 

where Oj{t) = ijq{nj{t)) includes the contribution due to the XPM effect. The expression 
for A2{t,z) can be derived by changing the indexes 1 <-> 2 in (34). The apphcation of the 
solution (34) requires the use of the extended algebra of time-dependent Bose operators that 
contains permutation relations of the form 

i,(ti)e^^(*2) = e^>(*2)+^^(*^-*i)i,(ti), (35) 
i.(ti)e^^(*^) = e^^(*^)+^(*^-*i)i,(ti), (36) 

where here T){t2 — h) — i^h{t2 — h) [22]. With expressions (35)-(36) for the average values 
of the Bose operators over the initial coherent summary state \cK{t)) we get: 

(e^^'W) = e*''^J'(*)"^^W, (e*^^'^*^) = e*'^^'^*^~'^^'^*\ (37) 
^eOj(*i)+Oj(*2)^ = ^^[<|yJitl)+<|yJ{t2)]-[^iJ{tl)+^^J{t2)]~lCj{tl,t2)^ ^3g-) 

^g(5,(ti)+(3j(t2)^ ^ e't^^'(*i^+^^(*2^]"t''^(*i^+''^'(*2)l"^^'(*i'*2). (39) 

The parameters 0j(t) = 2^nj{t), jlj{t) = 7^nj(t)/2 are connected with the XPM of the j-th 
pulse. Here (i)j{t) and ]Cj{ti,t2) are the nonlinear phase addition and the temporal correlator 
caused by XPM, respectively. 

The Stokes operators are defined in agree with the formulas (9)-(12) with D-matrices in 
the form: 

where z is the distance inside the anisotropic Kerr medium where the pulses overlap. By 
using Eqs. (35), (36) and Eq. (37) we estimate the average values of Stokes operators: 

{S2{t)) = 2[ni(t)n2(t)]'/'e-^i(*)-^^W cos [$2(t) - ^i(t)], (41) 
{S,{t)) = 2[ni(t)n2(t)]'/'e-^^(*)-^^W sin [$2(t) - $i(t)], (42) 

where here Aj{t) = [Xjif) + jijif), $j(t) = (pjif) — (pj(t) + ^j{t). In the approximation 
7j,7 ^ 1, by making use of Eqs. (35), (36) and Eqs. (38), (39), for the correlation function 
of 5*2 we get 
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i?52(^,^ + T) = 5(t) + /i(r)[rii(i)02(i) -ri2(t)(/)i(0]sin2[$i(i) -$2^] 

+ ^(r){ni(t)[0^(i)+0^(i)] +n2(i)[0?(i)+0?(t)]}sm^ - Mt)l (43) 

One can easily see that in the absence of XPM, i.e. 4>i{t) = 4>2{t) = 0, the expression above 
reduces to the expression (27) of the previous section. To derive the analytical expression 
of the fluctuation spectrum of S2 we apply the WK theorem. Thus, 

Ss,in,t) = l + 2L(fi)[ni(t)02(t)-n2(t)0i(t)]sin2[$i(t)-$2(t)]^ 

+ 4L\n){m{t)[<l)l{t) + 4>M + n2(t)[0Kt) + 0Kt)]}sin2 [$i(t) - ^t)]- (44) 

Notice that the XPM effect introduces specific contributions in the fluctuation spectrum 
[compare (44) with (29)]. Their change allows for an additional control of fluctuation spec- 
trum (see [26]). At the linear phase difference Aip(t) = (pi{t) — (p2{t) 

A^(t)„p. = larctan^ Mt)Mt) - Mt)Mt) 



2 \L{Qo){nom<i>m+m]+Mt)[m+m]}. 

+Mt)-Mt)-Mt)+Mt), (45) 

optimal at a deflnite frequency — ^^qTv, the spectrum (44) reaches the minimum 

s 8,(^0, t) = 1 + 2L\no) {MtMit) + 02 (t)] + n,{tM{t) + ^M) 

-2L{no) [[ni(t)02(t) - n2{t)Mt)? 

+L2(Qo)(ni(t)[0^(t) + 0^(^)1 +no,2(i)[0?(i) +0?(i)])']'^'. (46) 



In Fig. 5 we display the fluctuation spectra of 152 for different relations between the intensities 
of pulses. With the increase of intensity of the pulse 2 (the control pulse) a better level of 
suppression of quantum fluctuations can be achieved. A similar result, now displayed in 
Fig. 6, can be obtained by increasing the nonlinear coefficient 72 in comparison with 71. In 
conclusion, besides the optimal linear phase arrangement (45) that produce the suppression 
of quantum fluctuation at the frequency of interest, the increase of the intensity of the 
control pulse (Fig. 5), or of the nonlinear coefficient 72 (Fig. 6), can be employed for the 
achievement of desired level of squeezing. 



VII. INTERFERING QUADRATURE-SQUEEZED PULSES 

In this section we show that by overlapping a coherent pulse with an interference USP, 
obtained by mixing two quadrature-squeezed USPs on a BS, one can produce the simultane- 
ous squeezing in the first two Stokes parameters 5*0 and 5*1, as well as the squeezing in one of 
the last two, ^2 or ^3. We shall compare our theoretical results with the experimental ones 
presented in [5,6]. We start by analyzing the situation depicted in Fig. 7. Two quadrature- 
squeezed USPs fall simultaneously on an inclined BS with reflectivity R and transmissivity 
T, giving rise to reflected and transmitted interference USPs whose amplitudes depends on 
R and T. The outgoing nonclassical interference USP at the exit 1 of BS {Bi{t)) overlaps 
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a 

FIG. 5. Normalized spectral variance Sg^{^,t) at (f)o,i = 2 for initial phase difference Aip{t) 
chosen optimal at = 0. Curves are calculated at time t = 0, 71 = 72/4 = 2j, and correspond to 
"0,2 = no,i/4 (a), no,2 = no,i/2 (b), rlo,2 = "0,1 (c), no,2 = 3no,i (d) (from Ref. [26]). 




FIG. 6. Normalized spectral variance Sg^{Cl,t) at ^0,1 = 2 for initial phase difference A(p{t) 
chosen optimal at Qq = 0. Curves are calculated at time t = 0, no.i = ^0,2, 7 = 7i/2 and 
correspond to 72 = 271 (a), 72 = 871 (b), 72 = 471 (c), 72 = 871 (d), 72 = 671 (e), 72 = 771 (f) 
(from Ref. [26]). 
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FIG. 7. Interference on a BS of two independent quadrature-squeezed pulses produced by using 
the SPM effect in nonlinear electronic Kerr media (x^^^)- The outgoing nonclassical pulse at the 
output 1 (Bi(t)) is overlapped spatially with a coherent pulse field {A^lt)). 



with a third coherent USP {A^it)). We consider R and T being independent on frequency, 
direction of propagation or polarization and we treat both reflections symmetrically. In this 
latest case the relationship between input and output states on BS reads [27] : 





A2{t,l2] 



where Bj{t) is the annihilation operator on the output j of the BS and Aj{t, Ij 



(47) 



(j = 1,2). The Stokes operators at the output 1 are defined in agree with the formulas 
(9)-(12) when now the D-matrices are 



Bi{t) 
A-,{t) 



D2{t) 



B^{t) 



(48) 



By using the algebra of time-dependent Bose operators developed in [20,21], and summarily 
exposed in Section II, for the average values of Stokes parameters on the summary coherent 
state \(y.{t)) we find: 



(5o,i(t)) = i?ni(t)+Tn2(t)±n3(t) 

+ 2/RT[ni(t)n2(t)]^/^e-''i(*)-''2(*) sin [$2(t) - ^i{t)i 
{S2{t)) = 2VT[n2{t)ns{t)]'/'e-^^^'^ cos [^^it) - $2(t)] 



+ 2x/i?[ni(t)n3(t)]^/'e-'^iW sin [ip^it) - $i(t)]. 



(49) 



(50) 



where in (49) the sign '+' stands for {So{t)) and '— ' for (S'i(t)). The average value of 5*3 
can be easily obtained by shifting (50) in phase with — 7r/2. Notice the appearance of the 
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sin-term in (49) due particularly to the interference on BS. We calculate now the correlation 
functions of Stokes parameters in the approximation 7^ ^ 1 (j = 1,2). Below we write 
down the results obtained in the frame of the quantum model presented in Section II; 



i?5o.i(i, t + T)^ S{t) - h{T){2Vm'[ni{t)n2{t)nRMt) ± TMt)] cos [$i(i) - $2 W] 

+ RT[ni{t)cj)2{t)-n2{t)Mt)]sin2[^i{t)-^2{t)]} 
+ giT){RT[n,it)(j)lit) + n2(t)0?(t)] cos^ [$i(t) - $2(t)]}, (51) 
Rs,{t,t + T) ^S{T)+h{T)ns{t){RMt) sin2[<l>i(t)-(^3(t)]-T02(t) sin 2[<l>2(t) -<^3(t)]} 

+ g{T)m{t){Rcf>l{t) cos^ [$i(i) - <p3{t)] + Tcf>l{t) sin^ [^^{t) - <P3{t)]}. (52) 

The presence of the additional terms besides (5-function indicates that the formation of the 
PS state in all three Stokes parameters Sq, Si, §2 is permitted. The fluctuation spectra are 
straightforward with the use of WK theorem: 

Sso,A^,t) = 1 - 2L{Q){2^[ni{t)n2{t)]'/^[RMt) ± Thit)] cos [$i(t) - $2(t)] 

+ i?T[ni(i)02(i) - ^2(t)0i(t)] sin2[$i(t) - $2(i)]} 
+ AL''{n){RT[n,{t)(l>l{t) + n2{t)(j>i{t)] cos^ [^,{t) - $2(i)]}, (53) 
Ss,{^,t) ^l + 2L{n)ns{t){RMt) sm2[^,{t) - <ps{t)]-Tcf>2{t)sm2[^2{t) - <Ps{t)]} 

+ AL\n)ns(t){R(l)l(t) cos' [$i(i) - ^s(t)]+T(l>l(t) sin' [^2(t) - ^3(t)]}. (54) 

Since the correlation function of 5*3 and its fluctuation spectrum can be easily obtained by 
shifting (52) and (54) in phase with 7r/2, respectively, their analytical expression are not 
presented here. Notice that the above correlation functions and the fluctuation spectra are 
determined by the parameters R and T. In what follows we consider the most common case 
of a 50/50 BS, i.e., R = T= 1/2. 

We analyze now the first two Stokes parameters and by addressing the issue 
concerning the optimization of the linear phase difference Aip{t) = (pi{t) — 972 (^)- However, 
due to the analytical complexity of the expression (53), we restrict to the case ni{t)(f)2{t) — 
^2(^)0i(^) when sm-term in (53) vanishes. Thus, the resulting expression can be minimized 
if the linear phase difference Aip{t) at a defined frequency flo has the form 



A(^opt (t) ~ arccos 



Rni{t)±Tn2{t) / ni{t) 
2[ni(i)+n2(t)]0(t)L(Qo) \RTn2it) ^ 



1/2 



-Mt) + Mt)- (55) 



With this adjustment (55) for the linear phase difference between pulses the expressions (53) 
achieve the minimal value 

Sso,A^o,t) = 1 - [Rni{t)±Tn2{t)]'' [ni{t) + n2{t)]-\ (56) 

Thus, at any frequency ft the spectra (53) take the form 

<^ (Of^ <^ (o f^ I [RMt)±Tn2{t)]'[m-L{no)]^ 

Ss,Mt) = Ss,Ano,t) + — -^^y^-^^^ JJI^^ ■ (57) 

The fluctuation spectra of 5*0 and 5*1 for various relations between pulses intensities, at 
Aif{t) optimized at fig — 0, are displayed in Fig. 8 and 9, respectively. We conclude that 
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FIG. 8. Normalized spectral variance Sg^{i}, t) for linear phase difference A(/?(t) chosen optimal 
at = 0. Curves are calculated at time t = 0, R = T = 0.5 and correspond to n2,o = l-5ni^o (a)) 
n2,o = 2ni,o (b). 

the adjustment of the control pulse intensity can be employed to obtain a desired level of 
squeezing of quantum fluctuations in Sq and Si. This conclusion remain also valid in case 
the linear phase difference is optimized at frequency Qq = 1, the simultaneous suppression 
being correspondingly achieved at high frequencies Q ^ 1 {u ^ 1/Tr) [Figs. 10 and 11]. With 
the simplification introduced by the condition ni(t)02(^) = 'n2{t)4>i{t), while the suppression 
achieved in So is large, the one achieved in Si is small. 

However, it is important to remark that in the correlation functions and fluctuation 
spectra of 5*0 and 5*1 there is no information about the third coherent USP. As a consequence, 
the coherent USP can be disregarded in this respect or can be replaced with coherent vacuum 
which have an uncertainty identical to that of the coherent state. Basically, it means that 
the outgoing interference USP is already with squeezed quantum fluctuations in So and Si 
before the overlap with the coherent pulse. This situation due to the quantum interference is 
basically exploited experimentally in [5] where coherent vacuum is mixed on a polarization 
BS with a quadrature-squeezed USP obtained by using the SPM effect in an optical Kerr fiber 
(x^). Indeed, the outgoing interference USP shows the squeezed quantum fluctuations of 
about —3.7 and —3.6 dB in So and ^i, respectively. The simultaneous polarization squeezing 
of about —3.4 dB in both So and Si is also present when two quadrature-squeezed pulses 
interfere on a polarization BS [5,6]. A similar squeezing effect of about —3 dB in 5*0 and 5*1 
was reported in [8] where two bright quadrature-squeezed pulses obtained by employing the 
optical parametrical ampliflcation (x^) are mixed on a 50/50 BS. 

Let us now consider the next two Stokes parameters S'2 and S3. We focus on the param- 
eter S2 since once the suppression of quantum fluctuations is realized in it, the parameter S3 
becomes anti-squeezed, and viceversa. The fluctuation spectrum of Stokes parameter 5*2 can 
be easily minimized in the case 0i(t) = 02(t) = and — ^P2{t) = Tr/2. Actually, it 

was noticed experimentally in [9] that this relative phase of 7t/2 is critical for the detection 
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FIG. 9. Normalized spectral variance t) for linear phase difference Aip{t) chosen optimal 

at r^o = 0. Curves are calculated at time t = 0, R = T = O.b and correspond to n2,o = l-5ni,o (a)) 
n2,o = 2ni,o (b). 




FIG. 10. As in Fig. 8 but for = I. 
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FIG. 11. As in Fig. 9 but for = 1 

of squeezing. The linear phase difference A.(fi{t) — (fi{t) — (fi3{t) that at a defined frequency 
flo fulfills the condition 

cos2[0(t) + Aip{t)] = (i? - T)(f){t)L{Qo) sm2[(f){t) + A(p{t)] 

minimizes the expression (54). At this optimal linear phase difference the fluctuation spec- 
trum becomes 

Ss,{^0:t) = l + 2fi3(i)0'(i)L'(Jlo) -2ri3(t)0(t)L(fio)[l + (i?-r)VWi^'W]'/' (58) 
In consequence, at any phase Q the normalized spectra is 

Ss,ifl,t) = Ss,{no,t) + 2n3{t)(P\t)[L''{Q) - L\^^)\ + 1n^{t)m\^{^) - ^(^o)] 

X {1 + (i? - T)^^{t)L{Sl^){L{Sl) + L(Oo)]}[l + (i? - Tf<^\t)L\^^)\-^l\ (59) 

In Fig. 12 we displayed the fluctuation spectrum of S2 at time moment i = for different 
values of the nonhnear phase 0(i = 0) = 0o- The increase of 0o does not produce a better 
level of squeezing but moves the squeezing from low frequencies Q to high frequencies 
^"^^X. At the same time the squeezing at low frequencies is deteriorating. 

Note the obtained suppression of quantum fluctuations due to the stable spatial overlap 
of two quadrature-squeezed pulses of about of —2.8 dB in 5*2 reported in [5,6]. Indeed, since 
the parameter Si becomes squeezed, the parameter 5*3 is detected anti-squeezed with about 
-1-23.5 dB above the shot noise level. 

Summarizing, the mixing of two quadrature-squeezed pulses on a 50/50 BS followed 
by the spatial overlap of the outgoing interference pulse with a coherent field permits the 
simultaneous squeezing of quantum fluctuations in the first two Stokes parameters and 
(Si, as well as the squeezing in one of the remaining two, or S-^, (here S'2)- 
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FIG. 12. Normalized spectral variance Sg^{i^,t) for linear phase difference A(^(t) chosen opti- 
mal at r^o = 0. Curves are calculated at time t = 0, R = T = 0.5 and correspond to 0o = 0.5 (a), 
cPo = 0.75 (b), <Po = 1 (c), <po = 1.25 (d). 

VIII. CONCLUSION 

In this work we applied the quantum model developed in [22] to analyze the spatial 
overlap and interference of quadrature-squeezed USPs. We showed that the spatial overlap 
of coherent and quadrature-squeezed pulses or of two quadrature-squeezed USPs obtained 
by employing the SPM effect in the Kerr medium produces the squeezing in one of the last 
two Stokes parameters §2 or §3. We revealed that the adjustment of linear phase difference 
between pulses leads to the suppression of quantum fluctuations of 5*2 or 5*3 at the frequency 
of interest. Moreover, we found that the increase of the intensity of the control USP produces 
a better level of squeezing. 

By investigating the spatial overlap of nonclassical pulses in a nonlinear anisotropic Kerr 
medium in the presence of XPM effect we established that the increase of the intensity of 
the control pulse or of one nonlinear coefficient (72) in comparison with the another one (71) 
is able to realize the desired level of squeezing in iS'2 or S3. 

We studied the overlap of a coherent pulse with a nonclassical pulse produced by in- 
terfering two quadrature-squeezed USPs on a beam splitter. By studying the fluctuation 
spectra of Stokes parameters we showed that the nonclassical pulse exhibits the simultaneous 
squeezing of quantum fluctuations in Sq and Si. Besides, we found that its spatial overlap 
with a coherent pulse field produce the squeezing in 5*2 or S3. 
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